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Gluon radiation is an important mechanism for parton energy loss as the parton traverses the 
quark gluon plasma (QGP) medium. We studied the gluon emission in QGP using AMY formalism. 
In the present work, we obtained gluon emission amplitude F(h,p,k) function , which is a solution 
of the integral equations describing gluon radiation including Landau-Pomeranchuk-Migdal (LPM) 
effects, using iterations method. We define a new dynamical scale for gluon emission denoted by x. 
The gluon emission rate is obtained by integrating these amplitude function over h. We show that 
these obey a simple scaling in terms of this dynamical variable x. We define the gluon emission 
function g{x) for gluon radiation for the three processes g — >■ gg,q — >■ gq and g — s- qq. In terms of 
this g{x) function, the parton energy loss calculations, due to medium induced gluon radiation, may 
become simplified. 

PACS numbers: 12.38.Mh , 13.85. Qk , 25.75.-q , 24.85. +p 

Keywords: Quark-gluon plasma, gluon radiation, Landau-Pomeranchuk-Migdal effects, bremsstrahlung, 

gluon emission function, parton energy loss in QGP medium 






Oh 
I 

CD 



> 

vn 

o 
o 



X 



Lattice quantum chromo dynamics (LQCD) calcula- 
tions 'l| predict a transition from confined state in 
hadrons to a deconfined state of quarks and gluons above 
a temperature of 170 MeV or an energy density above 
IGeV/fm'^ . In the relativistic heavy ion collisions at 
RHIC at BNL with an energy density above 5GeV/fm'^, 
experimental measurements of several observables indi- 
cate such a transition to a deconfined state of matter 
[2-[al J for details see the reviews [a-Q • It is currently be- 
lieved that this deconfined state consists of a strongly in- 
teracting Quark Gluon Plasma (sQGP), behaving nearly 
like a perfect Uquid Q- Among these observables, jet- 
quenching phenomenon is so far an important signal for a 
hot dense medium formed. Jet suppression has its origin 
in parton energy loss in the quark matter by gluon radia- 
tion, which distinctly differs from energy loss in hadronic 
matter. For example the suppression of high-pr pions, 
from 3GeV to lOGeV, of BNL experiments can be ex- 
plained by assuming a deconfined state. 
We studied gluon radiation problem as this has direct 
application to the energy loss of partons while traversing 
the QGP medium due to the gluon bremsstrahlung pro- 
cesses. In addition to radiation, the elastic energy loss 
of partons traversing the QGP medium is important for 
heavy quark quenching, observed in RHIC experiments. 
For an exposition of theoretical and experimental results 
on parton energy loss the readers may see an excellent 
review [lO|. 

In the present work, we study the gluon radiation mech- 
anism. The coherent radiation processes involve multi- 
ple scatterings of the partons in the QGP medium dur- 
ing the gluon formation time. This leads to interfer- 
ence effects known as Landau-Pomeranchuk-Migdal effect 
(LPM). Gluon emission is discussed widely in literature 
[il|-[i4|- These works treated the parton energy loss on 
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the basis of avarage ener gy l oss dependeing on the path 
length. As emphasized in [18| , the bremsstrahlung (gluon 
emission), is characteristically different in the sense that 
it is a stochastic process. Starting with a group of par- 
tons of fixed energy, the bremsstrahlung process results 
in a broad spectrum of final partons of width compa- 
rable to its mean energy loss. Further, the LPM ef- 
fect has different parametric dependence on energy for 
soft and hard parts of emitted gluon spectrum. As com- 
pared to the case of bremsstrahlung photon emission, the 
gluon emission also involves an enhancement mechanism 
when the emitted gluon and quark are nearl y c ollinear, 
thereby a need to consider ladder diagrams (l6|. How- 
ever, unlike the emission of electromagnetic quanta, the 
emitted hard gluon feels the random colored background 
field. The resummation of these ladder diagrams leads 
to Swinger-Dyson type integral equations. In this work, 
we follow the formalism given in [l6[-[18| which imple- 
ments LPM effects by resumming the ladder diagrams. 
For calculating parton energy loss arising from gluon radi- 
ation, one needs the differential gluon emission rates and 
this is given by Eq[T] [18] in terms of the F(h,p,k) func- 
tion. The bremsstrahlung integral equations determine 
the gluon emission amplitude F(h,p,k) function given by 
Eql^l Here, the two dimensional vector h = p x k is 
of the order of 0{gT^). It is a measure of coUinearity 
and its magnitude is small compared to pk. The term 
SE{h,p, k) in Eql5]is the energy differential between ini- 
tial and final states. Here, m^ = m\,/2 = 2g^T'^ /?> and 
other quark thermal masses are m? = g^T^ /?>. This for- 
malism is very similar to the photon emission integrals, 
however, as mentioned before, the emitted gluon has color 
and therefore interacts with other scattering centers as 
well as soft background fields [l^. Accordingly, there 
are three terms in the integral equations of [5] involving 
collision kernel C(h). A typical ladder diagram for gluon 
emission is shown in Figure[TJ We solve this integral equa- 
tion for F(h,p,k) function, by using iterations method as 
discussed in [1^. The F(h,p,k) distributions for various 
values of parton and gluon momenta (p,k) were obtained 




FIG. 1: gluon radiation processes that contribute at order Os 



considering the mechanism g ~> gg, q ~> gq and g -^ qq 
processes using relevant factors {Ng — 2,ds = dyi(= 8), 
and Cs = Ca{= 3),Cf = 4/3). Fig 12] shows a few 
F(h,p, fc) distributions for these two processes (for high 
parton momenta) for various values of parton and gluon 
momenta (p,k). The real part is shown in figure (a) and 
negative of imaginary part in figure (b). Figl^Jc) shows 
the real part of the distributions for g — gg process for 
high and low values of incoming and outgoing gluon mo- 
menta (p,k). The calculations for g — >■ gg process are 
shown in FiglH In all these figures, p always stands for 
the incoming parton momentum. 
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I. GENERALIZED EMISSION FUNCTIONS FOR 
GLUON EMISSION 

Before we discuss the emission functions for gluon emis- 
sion, it is very instructive to recall the emission functions 



for photon emission. In our previous works [T^ [23] , we 
showed generalized photon emission function by integrat- 
ing (EqlS)) the corresponding p±-f(p±) distributions (see 
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/ arc in general functions of {pQ,qQ,Q^,T,as} and 
therefore, we defined the generalized emission functions 
(GEF) g in EqlHl which are functions of only x variables. 
These GEF (g) are obtained from corresponding / values 
by multiplying with c coefficient functions given in |19| . 
As an example, in Figure HI we show the results for 
GEF for bremsstrahlung (FigHlJa)). The solid curve 
in (a) is the empirical fit to this emission function, 
given by Ea llll Consequent to finding the emission 
functions like those given by EqUTl we expressed the 
imaginary part of retarded photon polarisation tensor 
for any po,qQ,Q^,T,as values by using EqfT2l Using 
this approach, we obtained the phenomenological fits 
to virtual photon emission rates from QGP for ladder 
processes with LPM effects [19|. We provided simple 
phenomenological formulae which are useful in model 
calculations for experimental dilepton yields. Following 
the procedure of generalized photon emission function, 
we now try to obtain the generalized gluon emission 
functions. For this, we integrate these distributions, 
that have been shown in Figs l2l3l over the variable h, 
as given in EqUHl This quantity / is strongly process 
dependent and is a function incoming and outgoing par- 
ton momenta, plasma temperature and strong coupling 
strength as denoted by the variables p, k, T and as- The 
integrated quantity / is plotted versus k in Figure [S] for 
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FIG. 2: (a) Shows the real part of h distributions of the 
h • 5RF(h) for gluons. Various curves are for various parton 
momenta (p) and gluon momenta (k) values as mentioned in 
figure labels as (p,k). The distributions are obtained using 
iterations method. 

(b) Shows the imaginary part of h distributions of the h ■ 
S5F(h) 

(c) Shows the real part of h distributions of the h • RF(h) 
for pure glue process. Various curves are for various parton 
momenta (p) and gluon momenta (k) values as mentioned in 
figure labels as (p,k). 



FIG. 3: (a) Shows the real part of h distributions of the 
h ■ KF(h) for g ^t qq process. Various curves are for var- 
ious gluon momenta (p) and quark momenta (k) values as 
mentioned in figure labels as (p,k). The distributions are ob- 
tained using iterations method. 

(b) Shows the real part of h distributions of the h • 5RF(h) for 
g — >■ gg process. Various curves are for various gluon momenta 
(p) and quark momenta (k) values as mentioned in figure la- 
bels as (p,k). 

(c) Shows the imaginary part of h distributions of the h-Q'F(h) 



the process g — > gg. Figure shows / plotted for different 
values of p labeled on the figure. As seen in figure, the / 
values are scattered. Therefore, we defined a variable z\ 
as given in Eall4l We show / values versus z\ variable 
in figure [SI As seen in figure (a), it exhibits a linear 
behavior on log-log plot, extending over nine orders of 
magnitude. This apparently gives an impression that z\ 
is a good dynamical scale for the process g ^ gg- In 
order to examine this, we plot I j z\ in figure (b). As seen 
in figure (b), the I j z\ values are scattered and exhibit 
no useful trends, showing that z\ is not a dynamical 
variable for this process. Therefore, we now define the 
dynamical variable x and a function g, for g -^ gg 
process as given in Eas ll5ll51 In the FigfTTl we show the 
function g versus x. As seen in figure, all / values for 
different parton momenta merge. We fit this data with 



an empirical curve together with parameters as given in 
Eqs ll7ll9l We denote this function in EqllT] as gluon 
emission function for the process g — >■ gg. 

We carried out this for the processes q -^ gq and g — > qq. 
The / values for these two processes also donot exhibit 
any trends as a function zi variable, however, x remains 
a good dynamical variable. We show these results for the 
process q ^ gq in FigIS] versus x. The curve in this fig- 
ure is given by empirical fit and parameters in Eas l20l2"51 
Therefore, for q ^ gq , we define gluon emission function 
g{x) as given in Eql^Hl We performed these calculations 
for the process g — > qq and these results are shown in 
FigM The curve in this figure is given by empirical fit 
and parameters in Eg 1241291 Therefore, for the process 
g — >■ gg , we define gluon emission function g{x) as given 
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FIG. 4: (a) The dimensionless emission function t;^ (a;) versus 
dynamical variable xt defined in Eq[8l Six cases of temper- 
ature and coupling constant values considered are mentioned 
in figure labels in different type symbols. The symbols repre- 
sent the integrated values of px distributions as a function of 
{po,1o,Q^,T,as} values. The solid curve is an empirical fit 
given by Ea llll 
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FIG. 6: The integral of gluon h distributions for the pro- 
cess g — >■ gg. The plot shows J -r^-^h ■ 5RF(h) versus the 
variable zi. The different parton momentum (gluon p) values 
are shown on the figure and temperature of plasma is taken 
T=1.0GeV. 
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FIG. 5: The integral of gluon h distributions shown in pre- 



vious figures for the process g - 
emitted gluon momentum (k). 
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/ (In^h ■ sRF(h) versus function 



Here incoming parton mo- 
menta (in the present case, gluon momenta p) has not yet 
been integrated. The different parton momentum (gluon p) 
values are shown on the figure and temperature of plasma is 
taken T=I.OGeV. 



in Eql22 

After obtaining the gluon emission function g{x) for these 
three processes, we can perform p integrations required 
in Eq[Tl i.e. integration in terms of dynamical variables 
X. This will give us differential gluon emission rates. In 
the jet-quenching studies, one needs to estimate the en- 
ergy loss of high energy partons while traversing the QGP 
medium. In this problem, the differential gluon emission 
rates estimated by integrating over p variable, will have 
to be coupled in order to determine the differential energy 
loss. These results were already shown by [lOl, 1 16l4lg | ■ 
In the following, we examine the integrand of differential 
gluon emission rate of the EqUJ i.e., without p integra- 
tion, and we denote this in short by gluon rate, which 
should not be confused with the gluon emission rates af- 
ter integrating over p. We obtained the / values by inte- 
grating the h distributions J ■7^^h-5RF(h). The splitting 

are as given in [1^ . In the FigfTUl we 

show the / values (red curves), splitting function (blue 
curves) and the gluon rate (black curves) at three differ- 
ent p values for the process g — > gg. We show similar 
results for the process q -^ gq in FigHT] and for the pro- 
cess g ^ qq in FiglT^ 
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FIG. 7: The integral of gluon h distributions for the process FIG. 9; The integral of gluon h distributions for the process 



The plot shows f ,^ ^y h ■ 5RF(h) versus the new g 



qq. 



This plot shows J ,^ w h ■ KF(h) versus the new 



variable x. As before, different incoming parton momentum 
(gluon p) values are shown on the figure and temperature of 
plasma is taken T=1.0GeV. 



variable x. The different incoming parton momentum (gluon 
p) values are shown on the figure and temperature of plasma 
is taken T^l.OGeV. 
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FIG. 



The integral of gluon h distributions for the process 



q — >■ gq. This plot shows J -^^h • 5RF(h) versus the new 
variable x. As before, incoming parton momentum (quark p) 
values are shown on the figure and temperature of plasma is 
taken T^l.OGeV. 
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FIG. 10: Black line versus gluon momentum (k) shows 

the integrand (without p integration) of the differential gluon 



for the process g 



gg. The 
gg given 



emission rate of Eq[l] — J^r- 

blue curve represents the splitting function for g 

by |rp^p_,_j.| of Eq[T]and the red curve represents the integral 

value J ,^ \ 2h ■ 5RFs(h;p, fc). The incoming parton (gluon 
in this case) momentum p (labeled as pi in figure) for this 
figure is fixed p = 10, 20, 50GeV as indicated on the figure. 
Temperature of plasma is T=1.0GeV. 



FIG. 11: Black line versus gluon momentum (k) shows 

the integrand of the differential gluon emission rate of Eq[T] 

for the process q — >■ gq. The curves are as in Figure 
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In conclusion, the gluon emission in quark gluon plasma 
including LPM effects has been studied at a fixed tem- 
peratures and strong coupling strength. We defined a 
new dynamical variable x for gluon emission. Further, 
we defined gluon emission functions (GEF) denoted by 
g{x) for the processes g — >■ gg, q — > gq and g — > qq. 
We have obtained empirical fits to these GEF and pro- 
vide the functional forms and parameters for all the three 
processes. We compared the differential gluon emission 
rates (without p-integration) for these three processes. 
In terms of the GEF, we may calculate the differential 
gluon emission rates for these processes. These empirical 
formulae will be useful in calculations of parton energy 
loss by medium induced gluon radiation. 
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FIG. 12: Black line versus quark momentum (k) shows the 
integrand of differential emission rate of Eq[l] - 



process g — >■ qq. The details are as in Figure [101 
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